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Abstract 



We address the problem of general dissipative regularization of the quasilinear trans- 
' port equation. We argue that the local behavior of solutions to the regularized equation 

near the point of gradient catastrophe for the transport equation is described by the 
[ logarithmic derivative of the Pearcey function, a statement generalizing the result of 

I A.M.Il'in [13]. We provide some analytic arguments supporting such conjecture and 

test it numerically. 

^ ' Introduction 

^ ■ In this article we address the problem of shock formation in a general dissipative regu- 

' larization 
^ ! 

O ■ ut + a{u)ux = e [b{u)uxx + c{u)ul] + \bi{u)uxxx + ci{u)UxxUx + di{u)u^^ + . . . (1) 

of the quasilinear transport equation 

; Ut + a{u)ux = Q, a'(n) 7^ 0, M,a;GM. (2) 

H ■ 

Here e is a small positive parameter, the coefficient h{u) does not vanish. In such a study we 
were inspired by the Universality Conjecture of [3j concerning the universal shape of disper- 
sive shock waves at the point of phase transition from regular to oscillatory behavior. This 
universal dispersive shock profile is described in terms of a particular solution of certain gen- 
eralization of the Painleve-I equation (importance of this particular solution in 2D quantum 
gravity and the theory of Korteweg-de Vries equation was also observed in [101 [D The 
universality conjecture for solutions to the Korteweg-de Vries equation with analytic initial 
data was proved in [2|. Further numerical evidences supporting the universahty conjecture 
of [3] can be found in [1] . Another starting point for the present research was the remarkable 
result by A.M.Il'in (see the book [13] and references therein) describing the asymptotics of 
the generic solution to the equation 

Ut + a{u)ux = euxx (3) 
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at the point of shock formation in terms of the logarithmic derivative of the so-called Pearcey 
integral (see below the precise formulation of the Il'in's asymptotic formula). In both dis- 
persive and dissipative cases the leading term of the asymptotic formula essentially depends, 
up to few constants, neither on the choice of a particular generic solution nor on the choice 
of a particular generic perturbation. 

Our main goal is to generalize the Il'in's universality result from the equations 
to the more general cas^ of equations of the form ([1]). In the present paper we present 
the conjectural form of such a generalization and describe results of numerical experiments 
supporting its validity. 

The paper is organized as follows. In the first section we explain simple arguments 
suggesting that, for sufficiently small e solutions to the perturbed equation ([T]) can be ap- 
proximated by solutions to the nonlinear transport equation (|2]) up to the time of gradient 
catastrophe of the latter. In order to save the space we omit the terms of order and higher 
from the formulae; their contribution to the asymptotic expansions will be of higher order 
anyway. We then proceed to the precise formulation of the dissipative universality conjec- 
ture (see Conjecture Elbelow) describing the leading term of the asymptotic expansion at the 
point of shock formation. We also give heuristic motivations of this main conjecture. In the 
last section we present results of numerical experiments supporting the main conjecture. To 
this end we begin with the standard Burgers equation in order to test the numerical codes 
based on the finite element analysis. Then we proceed to a particular case of generalized 
Burgers equation comparing the numerical solution with the asymptotic formula. 
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1. Critical behavior in the generalized Burgers equation 

Consider the following class of nonlinear PDEs depending on a small parameter e > 

ut + a{u)u^ = e [b{u)u^^ + c{u)ul] . (4) 

The coefficients a{u), b{u), c{u) are smooth functions, a'{u) ^ 0. The class of equations is 
invariant with respect to arbitrary changes of the dependent variable 

u^u = j{u), nu)j^O. 

Using such transformations one can reduce dl]) to one of the two normal forms 

ut + uu^ = e [b{u)uxx + c{u)ul] (5) 

or 

Ut + a{u)u^ = e b{u)u^^. (6) 

^It was shown in |8J that the Il'in formula also works for certain dissipative perturbations of the shallow 
water equations. 
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We will study solutions u = u{x, t; e) to the Cauchy problem 

u{x,0;€) = F{x) (7) 

with e-independent smooth initial data. In the particular case b{u) = 1, c{u) = one arrives 
at the generalized Burgers equation 

ut + a{u)ux = e Uxx (8) 

thoroughly studied by A.M.Il'in (see the book [13] and references therein). Let us briefly 
summarize the main results of |13j . 

For simplicity let us assume the initial data F{x) to be a monotone function on the 
entire real line x G M. The first issue is the comparison of the solution u{x,t]e) to the 
Cauchy problem (j4]), ([7]) with the solution v = v{x,t) to the inviscid equation obtained by 
setting e to with the same initial data 

vt + a{v)v^ = 

(9) 

v{x,0) = F{x). 

The two solutions asymptotically coincide on finite intervals of the x-axis for sufficiently 
small time 

t; £:) — f (x, t)| — 7- for £ — )■ 0+, xE[xi,X2\, 0<t<ti. 

However, the lifespan of the solution v{x,t) is finite, due to nonlinear steepening, if the 
function a{F{x)) is monotone decreasing on some interval of real axis. In this case the 
solution to the inviscid equation is defined only on the interval [0, to] where 



to = min ( — - — . (10) 
xm V [a{F{x))]J ^ ' 

Assuming the minimum in (fTOj) attained at an isolated point x = xq to be non-degenerate 
one arrives at a point of gradient catastrophe of the solution f(x,t), i.e., the limit 

lim f (xo, t) =: vq (11) 

t^to, t<to 

exists but the derivatives Vx{x,t) vt{x,t) blow up at the point (xo,to)- Thus the solution 
u{x,t;e) to the Cauchy problem (j4]), ([7]), if exists, cannot be approximated by the inviscid 
solution. For the equation (Q the right asymptotic formula was found in [13]. In the present 
section we will derive a suitable modification of this asymptotic formula and present some 
heuristic arguments justifying its validity. In the next section we will also give numerical 
evidences supporting our conjectures. 

Let us begin with recollecting some basics from the method of characteristics for solving 
the inviscid equation ([9]). For t <to the solution to the inviscid equation can be represented 
in the following implicit form 

x = a{u)t + f{u) (12) 
where the function f(u) is invers^ to the initial data v{x,0) 

f{v{x,0))=x. (13) 



^If the initial data is not a globally monotone function then the representation ([T2|) works on every interval 
of monotonicity. 
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Let (xo,to) be the point of gradient catastrophe of the solution. As above, denote 

Vo = v{xo,to) 

the value of the solution at the point of catastrophe. The triple {xo,to,vo) satisfies the 
following system of equations 

xq = aoto + fo 

= a',to + f, (14) 
= ao'to + /o • 

Here and below the following notations are used 

, (da{v)\ „ f d^a{v)\ 

ao = a{vo), fo = f{vo), = I — - — 1 , = I ^ 1 etc. 

\ / v=vo \ / v=vo 

In the subsequent considerations we will always assume that 

a', + 0. (15) 

The genericity assumption 

«::=-^K% + /o")7^0 (16) 

ensures that the graph of the solution f (x, to) has a non-degenerate inflection point at x = xq. 
Such a solution will be called generic. Locally a generic solution can be approximated by a 
cubic curve. For our subsequent considerations this well known statement can be presented 
in the following form (cf. [3]). 

Lemma 1 Near the 'point of gradient catastrophe a generic solution (fT2l) to the inviscid 
equation (Q admits the following representation 

v{x, t)=vo + k^/^v{x, t) + , A; ^ 0, t < to (17) 

X — xo — ao{t — to) T t — to , . 

^ = k ' ^ = ^ (18) 

where the function v{x,t) fort < is defined as the (unique) root of the cubic equation 

x = a'Qvi—KV^. (19) 

Proof can be easily obtained by substituting ( IT7|) . ( IT8|) into implicit equation (fT2l) of 
the method of characteristics and then expanding with respect to the small parameter k^^^. 
Observe that uniqueness of the root of the cubic equation (IT^ for t < is ensured by the 
condition 

a'oK > (20) 
valid due to a monotone decrease of the superposition a (f (x, to))- 



Remark 2 Observe that the cubic equation f[T9|) has a unique root also for t > to provided 
validity of the inequality 

m 2 /a'^^"^ 
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From this observation it is easy to derive existence and uniqueness of the solution v{x,t) to 
([9]) also for sufficiently small t — to > away from a cuspidal neighborhood 

\x - Xq - aojt - to)\ ^ ^ , ... . . n foo\ 
-, rr-p: < C tor some positive constant C [22) 

[t - tfiyy^ 

of the point of catastrophe. 

We are now in a position to formulate the main statement of the present paper. 

Conjecture 3 Let v{x,t) be the solution to the inviscid equation (Q with a smooth mono- 
tone initial data v{x,0) defined on M x [0,to) having a gradient catastrophe at the point 
{xojto) satisfying ( IT5l) and ( l20l) . Assume the smooth function b{u) to be such that 

bo := bivo) > 0. (23) 

Then 

1) for sufficiently small e > there exists a unique solution u{x,t;e) to the generalized 
Burgers equation with the same e-independent initial condition 

u{x,0;e) = v{x,0), xeR 

defined on M x [0,to + ^{^)) for some sufficiently small 6{e) > 0; 

2) away from a cuspidal neighborhood of the point of catastrophe the solution u{x,t;e) can 
be approximated by the inviscid solution v{x,t) 

\u{x, t; e) — f (x, t)\ = 0{e). 

For arbitrary X , T there exists the limit 

u (xo + ao/3 e^''^T + a e^'^, to + /3 e^/^T) - vo 
hm — ^ =: U(X,T) 24 

where 

•={?)"■ •=©)"■ -(A)'" 

Moreover, the limit does not depend on the choice of solution neither on the choice of the 
e-terms in the generalized Burgers equation (jl]). It is given by the logarithmic derivative of 
the Pearcey function 

U{X, T) = -2^ log r e-K-^-2-^^+4.x)^^_ (26) 

A somewhat stronger version of the last statement of the Main Conjecture can be given 
in the form of the following asymptotic formula 

/ , N , l/4rr — Xo — ao(t — to) t — to\ ,^ / l/2\ 

u{x,t;e)=vo + 7e ^ U \^ , j + [e ^ ) (27) 

expected to be true on some neighborhood of the catastrophe point. For the particular case 
b{u) = 1, c{u) = the asymptotic formula fl27|) coincides with the one obtained by A.M. Il'in 
(see in [13] )• 
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Let us add few heuristic motivations of the Main Conjecture. First, let us consider the 
small time behavior of the solution u{x,t;e). As the function v{x,t) satisfies (j4]) modulo 
terms of order e, one can seek the solution to the generalized Burgers equation in the form 
of a perturbative expansion 



u 



{x,t; e) = v{x,t) +ev'^'^\x,t) + e'^v^'^\x,t) + 



The terms of the expansion have to be determined from linear inhomogeneous equations (see 
details in [I3]). For example, the first correction can be found from the following PDE 

vi^^ + {a{v)v^^^)^ = b{v)v^^ + c{v)vl 

Instead, one can apply the method of the so-called quasitriviality transformations [5], [H] 
finding a universal substitution 

v^u = v + Y.''^'^'' ' ' .CT'^'^ ''^'^ (28) 

k>i "^^ 

transforming any monotone solution of the inviscid equation dH]) to a formal asymptotic solu- 
tion to the perturbed equation Here fk{v] v^-, v^x-, ■ ■ ■ , f log \vx\) are some polyno- 
mials in the variables v^, v^x, ■ ■ ■ , v^^''~'^\ log \vx\ with coefficients that are smooth functions 
of V. They satisfy the following homogeneity condition 

fk {v; X Vx, X'vxx, . . . , A^^-2^;(^'=~2), log \vx\) = X^'-^fkiv; Vx, Vxx, • • • , v^^'-'^ log\vx\), k>l 

(29) 

for any A 7^ 0. Advantage of the perturbative expansion written in the form (128|) is the 
locality principle: changing the unperturbed solution within a small neighborhood of a point 
{x*,t*) does not change the value of the perturbed solution outside the same neighborhood 
of the point. 

For convenience of the reader let us explain the computational algorithm for derivation 
of the perturbative expansion fl28|) . For simplicity let us consider a perturbed equation of 
the form 

Ut + uUx = e^{u]Ux,Uxx,---) (30) 

where $(«; Ux, Uxx, ■ ■ ■) is a smooth function of its variable polynomial in jets Ux, Uxx etc. 
We will rewrite (130|) as an equation for the function x{u,t) inverse to u{x,t): 

xt = u-exu^ (u;—,-^,...y (31) 

The clue is in the following statement (see [H]) describing the perturbative solution to fl5T]) . 
Lemma 4 Define the function \I/(m; Xu, Xuu, ■ ■ • ) by the formula 

'i/{u]Xu,Xuu, ■ ■ ■) = J (^u; ^ 

Then the function 



dxu- (32) 



x{u,t) = x'^^\u,t) - e x^^\u,t) (33) 
such that 

xf^=u, a:« = vl/(n;xr,4t---) (34) 
satisfies the perturbed equation fl?T]) modulo terms of orders'^. 
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Proof immediately follows from independence from t of higher -u-derivatives of x^^^ : 

7^ — 7:. = 7^ 7:. — = 0ml for m> 1. 

dt (9m™ Sm™ dt 

Inverting the series (133|) one arrives at the needed algorithm. 

Corollary 5 Let v = v{x,t) be a solution to the PDF 

vt + vv^ = 

satisfying 0. Then the function 

u = v + ev^'^ [v;—,--^,...] (35) 

satisfies the perturbed equation fl30|) modulo terms of order . 

For the particular case of the generalized Burgers equation (jl]) the first terms of the 
quasitriviality expansion read 



u = V — e 



b Vtt ca' — ba' 



H -p2 v^\og\v^\ 

a' Vt a' 



+ {e^) . (36) 



It would be interesting to rigorously justify that, for sufficiently small e the above mentioned 
algorithm produces the asymptotic expansion of an actual solution to the generalized Burgers 
equation. 

Let us now consider the solution to (jl]) in a neighborhood of the point of catastrophe. 
After a change of variables in the equation (jl]) 

X — Xq — a^it — to) = C'l^x 

t-to = e^/H (37) 
1/4- 

u — vq = e ' u 



one arrives at the equation 



Another substitution 



ut + a'oU Us = hoUsx + [e^'^) . (3^ 



x = aX, t = l3T, u = -fU (39) 
reduces the leading term of fl38l) to the standard form of the Burgers equation 

Ut + UUx = Uxx 
provided the constants a, /3, 7 satisfy the constraints 

a',^ = l, &o4 = l- (40) 
The Burgers equation can be solved by the Cole-Hopf substitution 

^7(X,r) = -2Aiogiy(X,T) 
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where W = W{X,T) solves the heat equation 

Wt = Wxx- 

The Pearcey function 

/oo 
-oo 

clearly satisfies the heat equation. Let us check that, using this function in the substitution 

u = -2^^W{X,T) 

one arrives at the correct asymptotic expression of the function u near the point of catas- 
trophe 

X = a'oui - Ku^ + (e^/^) (41) 
(cf. eq. f|T9|) above). Indeed, rescaling the integration variable 

C = e'/'z 

we rewrite the expression for u in the form 

-2a 7 e''/*— log e ^ dC 

^■^ J -oo 

where 

^(C; x,t) = l (C^ - 2C^^ + 4C . (42) 

For t < to the phase function has a unique minimum at the point Co = Coi^yt) determined 
by the cubic equation 

x-xo- ao{t - to) = ^(t - to) Co - "Co- (43) 
Applying the Laplace formula to the Pearcey integral 

oo 



U 



e dQ = — =e (1 + [e)) 



•'3Co^-^ 



/3 

and using the obvious formula 

dS (Co (a;, t);x,t) ^ Co(a^,t) 
dx 2a 

one arrives at the following expansion 

u = ^e-'l\,{l + iD{e)). 

Substituting into the cubic equation ( l43ll yields ( 1411) provided the constants a, /3, 7 satisfy 
one more constraint 

-,=^- (44) 
Together with the constraints f HO|) this gives f l25|) . 
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2. Solving numerically the generalized Burgers equa- 
tion. Comparison with the asymptotic formula. 

In order to test the numerical algorithms we will begin with the standard Burgers 
equation. First, let us consider the Cauchy problem for the inviscid equation 

Ut + uUx = (45) 

u{x,0) = F{x) 
At the point of catastrophe one has 

xo = ao + F{ao)to, to = - ] uq = F(ao), F"{ao) = (46) 



(cf. eqs. f lT^ above). For the particular choice of the initial data F(x) = 5- the point 

1 + X'^ 

of the catastrophe can be located as follows 

xo = V3, to = uo = -. (47) 

For t > to the solution to the Cauchy problem is close to a discontinuous one. Indeed, 
it is well known (see, e.g., [16]) that the limit at e — of a smooth solution to the Burgers 
equation 

Ut + uu^ = eu^^ (48) 

is described by a discontinuous function on the (x,t)-plane. The curve of discontinuity 
X = s{t) of the limiting function is called shock front (the solid line on Fig. [1]). We will be 
computing the numerical solution to the Cauchy problem in a neighborhood of the shock 
front and comparing it with the Il'in asymptotic formula. Let us explain the algorithm used 
for determination of the shock front. 

Fixing a point t = t* we will select an array of values {x*} in some neighborhood of the 
curve X = s(t). We will evaluate the function u = u* at the points (t*,x*) with the help of 
the Il'in asymptotic formula using Maple for computation of the Pearcey function. 




Figure 1: Shock front. 
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In order to determine the shock front (see [TSj) let us use the Rankin-Hugoniot condi- 
tions 

^ = -(F(ai) + F(a2)), (49) 

where ai(t) and 02 (^) are determined by the equations of characteristics 

x{t) = ai + F{ai)t, x{t) = a2 + F{a2)t. (50) 

Differentiating fl50l) in t and taking into account fj49|) one arrives at a system of differential 
equations for the functions x{t), ai{t), 02 (t) 

nx 1 

- = -(F(aO + i^(a2)) (51) 

c/ai lF(a2) + F(ai) 
^ ~ 2 1 + F'(ai)t 

(ia2 1 -F(ai) — F(a2) 
^ ~ 2 1 + F'(a2)t 
The initial data for these equations have the form 

ai(to) = oo, a2(to) = oq, x(to) = xq, (52) 



where xq, to, uq are determined by eqs. ( H6|) . 

If the solution to the Cauchy problem fl51l)- (!52l) can be written in an explicit analytic 
form then also the shock front can be computed explicitly. Otherwise the system (!5T!) - (l52l) 

can be solved numerically. Observe that at t = to one arrives at an ambiguity of the form -. 

It can be resolved with the help of asymptotic expansions of the functions x{t), ai{t), 02 (t) 
near the point t = Iq. If ai{t) < < a2{t), t > to then for the characteristics ai(t), 02 (t) we 
have 

aiW-o-(^(t-to)j , «^W-o+(^(t-to)j . (53) 

The expansion of x{t) near t = to has the form 

x = Xo + F{ao){t-to). (54) 

So, for solving the Cauchy problem fl5Tl) - fl52|) we will solve the system of differential equations 
( ISTl) where we put t = to + At. Here At is the time step. We use the asymptotic values (l53|) . 
(15^ as the initial data, i.e. 

x{to + At) =Xo + F{ao)At. 
In order to control the computation the following identity will be used (see e.g. [T^) 

r-ai 



^(Fi(ai) + F2(a2))(ai-a2) = [ ' F{ 

^ Jan 



a)da. 
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Finite element analysis 

For solving the standard Burgers equation fj48l) we will use the finite element method (see 
[6], [Z], |9], [H], [12]) realized in the package FreeFem++ [1^. Since this package, strictly 
speaking, is not designed for solving spatially one- dimensional problems one can reformulate 
the original problem as a 2D one considering solutions depending on one space variable only. 
Let us assume that the 2D domain has the rectangular form 

fi = {(x,y):0<x<L,,0<y<Ly}, 

of the size x Ly and ^ Ly. 

We will impose the no-fiux boundary conditions at y = 0, y = Ly but no specific values 
of u at X = 0, a; = Lj, assuming that boundary values of u are fixed at some fictitious 
boundary of a wider region 



du 
dn 



du 

= — 

y=0,Ly 



= 0. (55) 

x=0,Lx 



Here n is the exterior normal to the boundary dfl, d/dt = d/dt + ud/dx. 
In the numerical experiments we will use the following initial data 

"1-0 = IT^ (56) 

For the time approximation the semi-explicit Euler scheme will be used. To this end we 
multiply the equation by a test function 6 integrating the resulting expression over the 
domain Q 

II ^^^^ ^ ^^^^ ^ jjeutre,.,y 

n n 
or, taking into account the boundary conditions 

m+1 _ yiri \ f du 

e + uV+^e + dxdy = / e^ds (57) 

r J J on 

n on 

The problem fl57|) in the weak formulation along with the initial conditions f l56|) is solved 
by means of the FreeFem++ package. 

Comparison of the numerical solution with the Il'in's asymptotic formula near the shock 
front X = s{t) for t = 1.54, e = 0.01 is shown on Fig. |2l The solid line shows the numerical 
solution obtained by the finite element method while the dashed one corresponds to the 
asymptotic solution (!27|) . On the right hand part of the figure the region near the catastrophe 
point xq, to, uq (see (147|) ) is zoomed in. 
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-6 -4 -2 2 4 6 ^ 1.6 2 2.4 2.8 3.2 ^ 



Figure 2: Comparison of a numerical solution to Burgers equation with the Il'in's asymptotic 
formula for e = 0.01. 

On Fig. |3] the difference between the asymptotic solution ui given by ( 127|) and the 
numerical solution up is shown in the logarithmic scale. The evaluation of ui and Mi;' is done 
for 

t* = 1.54, X* = {1.74, 1.75, 1.76, 1.77, 1.78}, 
e = {0.0025, 0.005, 0.0075, 0.01, 0.025, 0.05, 0.075, 0.1}. 
The average slope is 0.5175 with the expected value 0.5. 

ln{ui - UF)i 




-6 -5.5 -5 -4.5 -4 -3.5 -3 -2.5 ln{e) 



Figure 3: Numerical estimate of the truncation error in the asymptotic formula. 

During the computation we control the total mass as function of time. With the bound- 
ary conditions under consideration the total mass is a conserved quantity. So the conser- 
vation of the total mass is a good test of the quality of numerical simulations. The results 
for e = {0.1,0.01,0.0025} are shown on Fig. |1 On the interval [0,1.8] with e = 0.1 the 
relative error is 0.0024, for e = 0.01 the relative error is 0.0006 while for e = 0.0025 it drops 
to 0.0006. 
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M(n) 
2.8 

2.78 

2.76 - 

2.74 - 
2.72 

2.7 



e = 0.1 

e = 0.01 

e = 0.0025 



0.4 0.8 1.2 1.6 2 t 



Figure 4: Testing conservation of the total mass. 



Generalized Burgers equation 

Let us now proceed to a particular example of the generalized Burgers equation (jl]) 

ut + uu^ = e{uux)x- (58) 

completed by the boundary conditions ( l55l) and the initial data ( 156|) . Like in the case 
of the standard Burgers equation ( l48l) the semi-explicit Euler scheme will be used for the 
time approximation. The variational reformulation of the problem along with the boundary 
conditions reads 



/ / { — - — e + vj^vj^+^e + evj^vj^+^e, \ dxdy = / e^ds 

n an 



(59) 



The numerical solution of the problem ( 158|1 . in the weak formulation, with the initial data 
(1561) will be computed with the help of the FreeFem++ package. 

Like above, let us compare the results of the numerical simulations with the predictions 
given by the asymptotic formula (|27ll . On the Fig . [5] the solid curve shows the numerical 
solution while the dashed one is the graph of the asymptotic solution (|2711 . On the right 
hand part of the figure a neighborhood of the point of catastrophe xq, to, uq is zoomed in. 
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-6 -4 -2 2 4 6 



Figure 5: Comparison of a numerical solution to the generalized Burgers equation with the 
asymptotic formula (I27|l for e = 0.01. 

On Fig. E] the difference between the asymptotic formula uj and the numerical solution 
Up is shown in the logarithmic scale, for the following values of {x,t,e) 

t* = 1.54, X* = {1.74, 1.75, 1.76, 1.77, 1.78}, 

e = {0.01,0.025,0.05,0.075,0.1}. 
One can observe the average slope of 0.5221 against the expected value 0.5. 



ln{ui — Up) 
-2.6 

-2.8 

-3 

-3.2 

-3.4 

-3.6 




x = 1.74 

y^i^ X = 1.75 

X = 1.76 

x = 1.77 

X = 1.78 



-4.5 



-4 



-3.5 



-2.5 ln{e) 



Figure 6: Numerical estimate of the truncation error in the asymptotic formula (127|) for 
solutions to the generalized Burgers equation. 



As above we used the conservation of the total mass valid for our particular case (!58l) of 
the generalized Burgers equation as a test of validity of the numerical scheme. The results 
are shown on Fig. [7] for the values e = {0.1,0.01}. On the interval [0, 1.8] for e = 0.1 the 
relative decay is 0.0025, for e = 0.01 it drops to 0.0044. 
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M{u) 




2.72 - 



2.7 I > 

0.4 0.8 1.2 1.6 2 * 

Figure 7: Testing conservation of the total mass for the numerical solution to the generalized 
Burgers equation (1551) . 
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